Quantum Computing with Atomic Josephson Junction Arrays 
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We present a quantum computing scheme with atomic Josephson junction arrays. The system 
consists of a small number of atoms with three internal states and trapped in a far-off resonant optical 
lattice. Raman lasers provide the "Josephson" tunneling, and the collision interaction between atoms 
represent the "capacitive" couplings between the modes. The qubit states are collective states of the 
atoms with opposite persistent currents. This system is closely analogous to the superconducting 
flux qubit. Single qubit quantum logic gates are performed by modulating the Raman couplings, 
while two-qubit gates result from a tunnel coupling between neighboring wells. Readout is achieved 
by tuning the Raman coupling adiabatically between the Josephson regime to the Rabi regime, 
followed by a detection of atoms in internal electronic states. Decoherence mechanisms are studied 
in detail promising a high ratio between the decoherence time and the gate operation time. 



I. INTRODUCTION 

Josephson effects originate from a tunneling of the par- 
ticles in the condensed modes between two superfhiids 
and reflect the phase difference of the macroscopic wave 
functions between the superfluids. Initially discovered in 
the superconductors [1,2], Josephson effects have been 
studied intensively in trapped atoms both theoretically 
and experimentally [3,4]. In the atomic case, Joseph- 
son junctions can be constructed between two superflu- 
ids spatially separated by a double well potential and 
can be constructed between atomic internal modes cou- 
pled by lasers. Studies include the macroscopic quantum 
coherence between two atomic condensates and the ob- 
servation of the Josephson dynamics. [5] 

One important application of the Josephson junctions 
discussed in recent years is in quantum computing. Var- 
ious superconducting Josephson devices have been pro- 
posed for implementing a quantum computer, including 
the charge qubit, the flux qubit, the phase qubit. These 
qubits have been experimentally tested and have shown 
quantum coherent oscillations between macroscopically 
distinguishable states [6-8]. 

The atomic Josephson junctions can also be explored 
for quantum computing. In this paper, we present a 
candidate for implementing an atomic "flux" qubit with 
small number of atoms in an optical trap. We assume 
that a Bose Einstein condensate with three atomic states 
is stored in the lowest vibrational state of an optical 
trap [9]. The three internal atomic states correspond 
to three bosonic modes. Each mode is the analogue of a 
superconducting metallic island. Raman lasers generate 
the Josephson links between the internal modes, while 
atomic collisions provide an effective capacitive couplings 
between the modes. The phase differences between lasers 
plays the role of the magnetic field in the superconduct- 
ing loop. With the competition between the Josephson 
energy and the collision energy, the atoms behave collec- 
tively and the stationary states of the qubit have a co- 
herent particle transfer -the persistent current - between 
the internal modes. With only 15 atoms [10], the atomic 



counterpart of the superconducting flux qubit [7,8] can 
be realized, which bears all the qualitative features of the 
superconducting flux qubit. 

Compared with the superconducting flux qubit, the 
parameters of the atomic "flux" qubit can be controlled 
with large flexibility and high uniformity. Both the 
Josephson coupling and the collision interaction can 
be adjusted by external electro-magnetic sources. The 
Josephson couplings of different junctions can be made 
to high accuracy with the fine control of laser. While 
for superconductors, not only that the junction param- 
eters fluctuate due to the inaccuracy in fabrication, but 
the parameters are fixed for one sample. This advantage 
makes it easier to scale up the number of qubits in the 
atomic systems and provides various ways to implement 
gate operations. Another merit of the atomic qubit is 
that a projective measurement can be performed by adi- 
abatically switching the Raman couplings. On the con- 
trary, an efficient readout for the solid-state qubits is a 
problem many people are studying. The drawback of the 
atomic qubit is the slow gate speed which is limited by 
the strength of the collision interaction. Meanwhile, this 
drawback can be compensated by the long decoherence 
time. In the solid-state systems, various elementary exci- 
tations can damage the coherence of the quantum states 
in a time that is only one order longer than the gate time; 
while we show that in the atomic qubit, the decoherence 
time is thousand times of the gate operation time. 

In the following, the major results are summarized. 
In section II, we briefly review the superconducting flux 
qubit and the experimental achievement for the flux 
qubit. In section III, we give a detailed description of our 
proposal for the atomic "flux" qubit and how it can be 
realized experimentally. We also characterize the qubit 
at different parameter regimes and present typical en- 
ergy scales for the qubit. In section IV, we introduce 
the phase mode to compare this qubit with the super- 
conducting one and show that a small number of atoms 
indeed represents the macroscopic behavior of a Joseph- 
son junction. This section is followed by section V where 
the implementation of quantum logic gates are studied. 
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In section VI, a quantum nondcmolition measurement 
scheme is constructed via the adiabatic switching of the 
Josephson couplings. The decoherence of the qubit is dis- 
cussed in section VII. The conclusions are given in section 
VIII. 



II. THE SUPERCONDUCTING FLUX QUBIT 

Josephson junctions have been proved to be a promis- 
ing building block for quantum computers. Various 
proposals of Josephson circuits at different parameter 
regimes have been studied [6-8]. Among these, the super- 
conducting flux qubit — also named the persistent-current 
qubit — has been intensively studied both theoretically 
and experimentally. In the following, we briefly sum- 
marize the basic facts of the flux qubit in superconduct- 
ing Josephson Junctions, to allow comparison with the 
atomic "flux" qubit introduced in the following section. 

A. The Circuit of the Qubit 

The superconducting flux qubit [7,8] is a superconduct- 
ing loop with three Josephson junctions in series, as in 
Fig. la. Written in terms of the phase differences cross 
the bottom two junctions, ipi and <p 2 , the Hamiltonian 
(Eq. (11) in [8]) is 

Ht = ^P T M~ 1 P + Ej{2 + a - cos <pi - cos (p 2 

-a cos(2tt/i +<fi!- f 2 )} (1) 

where P = (Pi,A) T are the conjugates of the phase 
variables and have the physical meaning of the charges 
on the islands. The first term is the capacitive energy 
with M = (h/2e) 2 C where C is the capacitance matrix 
of the circuit. The rest of the terms are the Josephson en- 
ergy with Ej = I c {-^) 2 and I c being the critical current 
of the junctions. The third junction at the top of the cir- 
cuit has a Josephson energy of aEj with a = 0.75. The 
magnetic flux in the loop, f\ in unit of the flux quantum 
$o = h/2e, is an important control parameter for the 
qubit. Both the stationary states and the one-bit logic 
gates are controlled via this flux. 

The Hamiltonian in Eq. (1) describes a phase parti- 
cle in a two-dimensional periodic potential as is shown 
in Fig. lb. Each unit cell has two energy minima and 
is a double well potential. At f\ = 0.495, the lowest 
two states of the qubit localize in one of the two wells 
respectively and have opposite circulating currents. At 
fi = 1/2, the lowest two states are symmetric and anti- 
symmetric superpositions of the localized flux states and 
the energy splitting is due to the tunneling of the flux 
states over the potential barrier. Considering only the 
lowest states, the qubit can be described by the Pauli 
matrices for a 1/2 spin: H q — ^a z + ^o x , where the 
eigenstates of a z are the localized flux states and eo varies 



linearly with (/i — 1/2). Typically, the Josephson energy 
is Ej = 200 GHz and Ej/E c = 80. Numerical calcu- 
lations of the energy and current are shown in Fig. lc. 
The energy difference of the qubit states at f\ = 0.495 
is uo q <~ 10 GHz with the average currents of ±0.7J C ; at 
h = 1/2, to = 10 GHz. 

For a quantum circuit to be a good qubit for fault- 
tolerant quantum computing, five requirements have to 
be met [12]: 1. to identify a scalable quantum system; 
2. to perform universal quantum logic gates; 3. to pre- 
pare the initial state; 4. to read out the qubit states; 
5. to have a decoherence time longer than 10 4 quantum 
operations. The three-junction loop behaves as an effec- 
tive two-level system and can be mapped onto a 1 /2-spin. 
The qubit can be prepared to the ground state by cooling 
it to a temperature of T <~ 50 mK <C uj q . 

B. Quantum Logic Gates 

To achieve universal quantum logic operations, two el- 
ementary gates are required: single qubit rotation and 
two qubit controlled gate. For the superconducting flux 
qubit, the single qubit gate is by applying microwave 
oscillations to the superconducting loop. Typically, the 
Rabi frequency is uj r = 10 — 100 MHz in proportional to 
the amplitude of the microwave. The two qubit gate is 
constructed via the coupling of the circulating currents 
of the two qubits: Hint = Mi2\(Ii)(l2)\ with Ii, 2 being 
the currents of the two qubits and M\ 2 being the mutual 
inductance. The interaction can be of order of 1 GHz. 



C. Qubit State Readout 

The qubit is measured by inductively coupling the 
qubit to a dc SQUID magnetometer which is a supercon- 
ducting loop with two Josephson junctions as is shown 
in Fig. Id. When the current that flows through the 
SQUID increases, the SQUID stays in the superconduct- 
ing state until a critical current 1%^ where the SQUID 
makes a transition to a finite voltage state. The criti- 
cal current is varied by the flux generated by the qubit: 
51^? = ±S(p q I^ q simrf ex where f ex is the external flux 
in the SQUID and ±5ip q are flux of the two qubit states 
respectively. By measuring the critical current, the qubit 
states are read out. Due to fluctuations, the measured 
critical current has a distribution that is wider than 
Sl^ff which results in a nonprojective measurement of 
the qubit. 



D. Decoherence 

Many factors can result in quantum errors against the 
superconducting qubits. First, the errors can come from 
the imperfect control of the qubit circuits, for example, 
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off-resonant transitions during gate operations and un- 
wanted dipolar couplings between qubits. These errors 
can be prevented by the quantum control approach. Sec- 
ond, the fluctuations of the environment of qubit can 
cause decoherence of the qubit. In the solid-state qubits, 
many elementary excitations exist that can damage the 
qubit state, including the dipolar interactions between 
the qubit and the nuclear spins, the background charge 
fluctuations, and the noise coupled to the qubit from the 
measurement circuits. The decoherence time measured 
in experiments is lOOnsec [11] which is about 10 times 
of the operation time. This gives a lower bound for the 
generic decoherence of the qubit. 

III. THE ATOMIC "FLUX" QUBIT 

In this section we present an atomic counterpart of 
the superconducting flux qubit. The qubit is made of a 
mesoscopic Bose Einstein condensate of three- level atoms 
trapped in the lowest motional states of an optical trap, 
and interacting with each other via cold collision. Joseph- 
son junctions, which are the building block of this qubit, 
are constructed by laser coupling of the three bosonic 
modes of the trapped atoms. 

A. The Physical System and the Hamiltonian 

We consider a small number of three-level atoms 
trapped in a ID optical lattice, as is shown in Fig. 2a, 
which is described by the Hamiltonian 

Ho = ]T J dx^ a (x) (-^ + V(x)\ i> a {x) 

+ dx (n af3 {x)i)i(x)ip(}(x) + h.c) , (2) 

where the three internal states are labeled by a, (3. The 
first term is the single particle energy in a harmonic trap- 
ping potential: V{x) — ^moj 2 x 2 + ^mui 2 _(y 2 + z 2 ), where 
u; j^n are the trapping frequencies in the transversal di- 
rection and the longitudinal direction respectively. With 
a cigar-shaped geometry, we have uj± 3> 0J\\. We as- 
sume that the trapping frequencies are much larger than 
all other relevant time scales (e.g. the qubit energy and 
the gate speed) so that the atoms stay in the motional 
ground states and the qubit can be described by a three- 
mode Hamiltonian. The second term in Eq. (2) is the 
collision interaction. In practice, this term can be more 
complicated due to the scattering between the atoms in 
different internal states: for example, in the case of a hy- 
perhnc F = 1, Mp = 0, ±1 there will be collisional terms 
where two atoms in state Mp = collide to produce an 



atom in states Mp = +1 and —1 states. These terms 
can be suppressed by shifting the atomic levels by exter- 
nal fields, so that energy conservation in the collision is 
violated, i.e. these terms average away in the Hamilto- 
nian. Furthermore, detailed numerical studies show (see 
below) that we can simplify the collisional term to a fully 
symmetric interaction with U a /3a'/3' — 5af35 a a'5pf3'Uo, the 
main physical properties of the qubit are well preserved. 
Here, the interaction strength is 

m J 

where <f>(x) is the motional ground state of the trapping 
potential and a s is the s-wave scattering length with m 
being the mass of the atoms. With fixed number of 
atoms, the interaction strength increases with the density 
of the atoms. The last term in Eq. (2) is the Josephson 
couplings between the internal states generated by Ra- 
man transitions as is labeled in Fig. 2b as Q a /3- Both 
the amplitudes and the phases of these couplings can be 
accurately controlled by adjusting the laser parameters. 
In this system, we let Q ac = Qb c = and ^ a 6 = ^ie^° , 
where fii/f2 ranges between 0.5 and 1.5 and is an impor- 
tant factor for the speed of gate operations. The phase 
4>q is the analogue of the magnetic flux fi of the super- 
conducting qubit and is an effective controlling knob for 
the quantum logic gates. 

With the above discussion, the qubit can be well de- 
scribed by a three-mode Hamiltonian, 

n = U J2^+ J2 ( n ctpalaf, + n^ajaa) (4) 

a,/3 (a,f3) 

where N a is the number operator for the mode a. Here, 
the collision energy is the slowest energy scale which lim- 
its the speed of quantum logic gates, while the Raman 
couplings can be well controlled by lasers. In practice, 
Fcshbach resonances can be exploited to adjust the scat- 
tering length by several orders of magnitude [14,15] and 
the gate speed can be improved. 

The basis element in this qubit is to construct atomic 
Josephson junctions with small number of atoms. Atomic 
Josephson junction has three distinct parameter regimes 
[3]: (1) the Fock regime with U > ^oN t ; (2) the Joseph- 
son regime with U N 2 > £l N t > U ; and (3) the Rabi 
regime with floN t ^ UqN 2 . In the Fork regime, the col- 
lision energy dominates over the Josephson coupling and 
the eigenstates have fixed number of atoms in each in- 
ternal state. In the Josephson regime, the qubit behaves 
as a phase particle in the Josephson potential energy. 
In the Rabi regime, the atoms behave as noninteracting 
particles described only by the Josephson couplings. In 
a superconducting Josephson junction, the Rabi regime 
can never be approached with the enormous number of 
Cooper pairs. While for the atomic Josephson junctions, 
all three regimes are possible. In this paper, we assume 
the atomic qubit to be in the Josephson regime. When 
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compared with a large ensemble of atoms (say 10 5 atoms) 
in a superfluid state where the three-mode approxima- 
tion becomes inaccurate during fast gate operation, this 
system has the advantage that the three-mode model is 
robust against the qubit dynamics. 

In the Josephson regime, with N t ^> 1, Eq. (4) can be 
approximated by a phase model [3]. We introduce the 
phase variable (p a ,b that are the conjugate operators of 
the number operators N a ,b respectively. Due to particle 
number conservation, N c is not an independent operator 
with N c = N t — N a — Nb- By neglecting terms of order 
l/^/Wt, the Hamiltonian is 

T-Lphase = - ^2 U aPj. p. 1" -N t Q {cOS ip a + COS if b ) 

2 

+ -N t fli COs(f a - Ifib + 4>o), (5) 

which shows that with large number of atoms in the 
qubit, the Hamiltonian in the phase model maps ex- 
actly to Eq. (1) of the superconducting flux qubit with 
E c = 3C/ /4 [16], Ej = 2tt N t /3. In the next section, we 
will discuss the validity of the phase model with small 
number of atoms. 

We illustrate our model with the following parameters 
for 15 sodium atoms in the trap. For 23 Na, we choose 
the trap size to be L|| = 0.85 /im and L\\ = 10Lj_, which 
can be achieved with a far red detuned laser. The trap- 
ping frequencies are w\ = 3.7 kHz and oJq = 370 kHz. 
Let a s = 10 nm. With a density of p = 3 x 10 14 cm~ 3 , 
the collision interaction is Uq = 550 Hz. The Josephson 
couplings can be controlled so that ^l a /3\/N a Np >• Uo 
in analogy to the superconducting flux qubit. We let 
2il (N a ) w 70Uo in the following calculations. In the 
notations of the superconducting qubit, Ej/E c w 90. 



B. Effective Two Level System 

We have numerically studied the Hamiltonian in 
Eq. (4) with the above parameters. In Fig. 3a, we plot 
the energies and the average currents of the eigenstates 
of the qubit versus the phase (fio in the range 0.48 to 
0.52 (in unit of 2tt). The current operator is defined 
as I ac = i^o{a a a c — Sja a ). We define the lowest two 
states as the effective two level system of a qubit. At 
4>o = 0.495, the qubit energy is uo q = 1.4 kHz, where the 
states are labeled by the arrows. The stationary states 
have a coherent transfer of the atoms between the inter- 
nal states that provides a persistent particle current for 
the qubit. The currents of the two qubit states flow in 
opposite directions just as in the superconducting qubit, 
with (I ac ) i,2 = ±4.3£!o- This shows that the atoms be- 
have collectively just as the electrons in the supercon- 
ducting wires, which is a result of the interaction be- 
tween the atoms. For comparison, the energies of the 



qubit when Uq = are also plotted as the dashed lines 
in Fig. 3a. 

At fa = 1/2, the energy splitting to is 660 Hz which 
is the counterpart quantum tunneling in the flux qubit 
and an important feature of the qubit that is crucial for 
gate operations. We studied this splitting with various 
circuit parameters. Our result shows a dramatic depen- 
dence of to on the ratio between the Josephson couplings 
ro = fii/fio: at r o = 0.75, to = 1.1 kHz; at r = 1, 
t = 0.2 Hz; and for ro > 1, to is almost unchanged as r 
increases. 

In this study, we choose the fully symmetric interac- 
tion in Eq. (4) for the comparison with the supercon- 
ducting qubit. It can be shown that the detailed form 
of the interaction doesn't change the main feature of the 
qubit. For example, let the interaction take on the form 
V = C2F1 • F2, where C2 = (52— <7o)/3 and gi = 4irfr a^jM 
is the scattering strength in the F t = i channel, with Fi 
being the angular momentum of the atom i and F t be- 
ing the total spin [13]. Numerical results with C2 < 
shows that the calculated energy spectrum, and hence 
the properties of the qubit, has the same butterfly shape 
as that in the fully symmetric interaction, as is plotted 
as the dotted lines in Fig. 3(a). In this plot, we choose 
the interaction to be — Uq with CIi/Qq = 0.85. 



C. Finite "Size" Effect 

To see the effect of the small number of the atoms, we 
calculate the energies with various numbers of atoms, as 
is shown in Fig. 3b. The plot shows that the energies 
of the qubit converge as the number of atoms increases. 
Furthermore, it shows that when Nt = 15 the states of 
the qubit well represents the key features of a supercon- 
ducting flux qubit — the features of a qubit in the phase 
model. The surprising fact is that with a small num- 
ber of atoms, the atomic qubit reflects the properties of 
the flux qubit with over 10 10 Cooper pairs: the qubit 
states have opposite persistent currents; the phase in the 
Raman coupling induces energy difference that is nearly 
linear with (fi — 1/2; besides, even the wave functions in 
the phase space can be described by the localized phase 
states. 

The wave function in the basis of the phase vari- 
ables is \ifi) = J dtp a d<Pb\<Pa,tPb}( l Pa,tPb\i>}- In our cal- 
culation, we use the number state basis for the states: 

^n a ,rib 

\n a , rib), where c na . nb is the coefficient 
of the wave function. The wave function in the phase ba- 
sis is then: (9«,^W = E„ t ,„ l c,..,» l e" iw ""*" 1 ' In 
Fig. 4, \{<p a , </3h|V')| 2 °f the ground state is plotted in the 
phase basis with N t = 15, 30, 60 respectively. 

The phase model predicts that at (fio = 1/2 the wave 
function be a superposition of two local flux states. For 
the small number of atoms with a weak interaction, 
Fig. 4a shows that the qubit state localizes at the center 
of the phase space in contrast to the phase model pre- 



4 



diction; while, with N t = 60, the state is a superposition 
of two local states in agreement with that of the phase 
model. Fig. 4c shows the same result for fa — 0.495. 
With a stronger interaction, Fig. 4b and Fig. 4d show 
that the state of N t — 15 atoms agrees with the phase 
model result. Our study indicates the behavior of the 
qubit depends strongly on the factor U$Nf /QoN t . When 
UoNj- < SloN t , the qubit enters to the Rabi regime and 
single atom behavior starts to dominate over the collec- 
tive behavior. Our result also shows that with N t = 15, 
the qubit represents the main features of a phase-model 
qubit. 



IV. GATE OPERATIONS 

Now we are going to discuss how to realize quantum 
logic gates, qubit initialization, qubit state readout and 
the decoherence of the atomic flux qubit in the following. 



A. One-bit Gate 

The superconducting qubit is operated with external 
magnetic fields where microwave pulse in resonance with 
the qubit frequency is radiated on the superconducting 
loop. Off-resonant transitions to other states of the qubit 
can be neglected since the Rabi frequency is much smaller 
than the detuning. 

A similar scheme can be applied in the case of the 
atomic flux qubit. If we take a Raman laser coupling any 
two of our bare atomic states which make up the qubit, 
and we tune these lasers to match the energy difference 
of the qubit states, we can perform Rabi rotations be- 
tween the states. In order not to excite any higher lying 
states, the Rabi frequency should be less than the level 
spacings. In the atomic flux qubit, the qubit frequency 
and the detuning are of the order of 1 kHz which makes 
these gates slow. A first way to improve this, is to use 
adiabatic passage, i.e. a sweep of the detuning across the 
resonance, which allows a single qubit rotation on the or- 
der of the level spacing. Below we discuss in more detail 
another scheme based on fast switching of the phase fa 
of the Raman coupling f2 f,. 

Assume Ha = 7~io(fa = 0.495) and Hb = Ho(fa = 
0.5), and [Ha, Hb] ^ 0. We know from group theory 
that by switching the phase alternatively between these 
two phase values, any desired unitary transformation can 
be constructed within reasonable number of switchings as 
U = e- iHAt2 ™e- iHBt2 ™- 1 ■ ■ . e ~ iH At2 e ~iH B t 1 by a dj us ting 
the durations U of the pulses [17] . For a single qubit gate, 
we want the unitary transformation to be block-diagonal 
between the two qubit states and the other states. A nu- 
merical optimization of the {ti} is applied to a 12-pulse 
sequence of the Ha and Hb operators for the lowest six 
states of the qubit. We construct a NOT gate and a 



Hadamard gate Uh- The elements of the unitary opera- 
tors \Uij\ are shown in Fig. 5a, Fig. 5b. The off-diagonal 
elements t/»,i, f7»,2 <C 0.01 shows a high fidelity for one- 
bit gates. The total time for the gates is t\ ~ 2 msec for 
both gates. The accuracy of the gate can be improved 
by increasing the number of pulses in the sequence while 
keeping the total gate time short (which means faster 
switching of the operators Ha,b)- 



B. Two-bit Gate 



Two-bit gates can be constructed by external Joseph- 
son tunneling between neighboring qubits. As we men- 
tioned earlier, external Josephson tunneling is the tun- 
neling of atoms between spatially separated condensates. 
With the geometry in Fig. 2a where the qubits are aligned 
parallel along the longitudinal direction of the cigar- 
shaped trap, atoms with the same internal mode can 
tunnel from one lattice site to its neighbor site by switch- 
ing on a laser pulse for a short time. The tunneling is 
enhanced by a factor of N t of the number of atoms. 

We consider the tunneling interaction, 

H 2 = ttt^2 ( a ia fl 2" + a L«ia) ( 6 ) 

a 

the index 1 and 2 in the operators refer to qubits 1 and 2. 
The tunneling matrix can be estimated with WKB ap- 
proximation: fl t ~ tj^ exp (— AU/huj±), with u>±_ being 
the plasma frequency of the atoms in the trapping po- 
tential, and AU the trapping barrier for the qubit. The 
single particle tunneling Q t is enhanced by the number 
of particles and so does the speed of two-bit logic gates. 
The tunneling rate can be controlled by adjusting the 
laser pulse. 

The interaction H.2 can be calculated numerically. The 
matrix elements of the operator (ffljjy is obtained by cal- 
culating the overlap between the states |ijv t +i) f° r N t + 1 
atoms and the states o) a \jN t )- Our calculation shows that 
this interaction as well as that of the single-qubit gate in- 
duces coupling to the higher states of the qubits. This 
problem can be prevented by the same approach as that 
of the single-qubit gate — fast pulse sequence to decou- 
ple the lower states from the higher states. We apply 
a pulse sequence of 36-pulses with Ha = H.2(fa) and 

H B = H<£\fa) + H ( o\fa) where W< 1,2) are single qubit 
Hamiltonian at fa. In Fig. 5c, we show the absolute 
values of the matrix elements for the two-bit transforma- 
tions at fa = 0.495 and fa — 0.5 respectively. With a 
total pulse duration of 2 msec, the fidelity of the gates for 
N t — 15 atoms is higher than 99%. 
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V. ADIABATIC PROCESS AND 
MEASUREMENT 



B. Initial State Preparation 



The qubit we studied in the previous sections works in 
the Josephson regime where UqN^ Qc,N t ^> Uq. In 
this section, we present a quantum non-demolition mea- 
surement scheme during which the qubit is switched adi- 
abatically between the Josephson regime and the Rabi 
regime where £loN t ^> UaN^. In contrast to the mea- 
surement of solid-state qubit where it takes efforts to 
build efficient measurement schemes, our method pro- 
vides an easy-to- realize way for qubit readout. The same 
approach can also be applied to initialize the qubit. 



A. Qubit in the Rabi Regime 

In the Rabi regime, when the Josephson energy is much 
larger than the collision energy, we neglect the collision 
term and the qubit is described by the single atom Hamil- 
tonian, 







Hj = ( at, b\ ct ) \ Q 





n \ 
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which describes a three-mode atom where the internal 
modes are coupled by lasers. The eigenstates can be de- 
scribed by atomic states as 



(8) 



i=i 



where Sj and Si are the operators for the atomic eigen- 
states and ei are the eigenenergies with t\ < e 2 < and 
2(f2 — £i) < (^3 — The ground state and lowest ex- 
cited states of the qubit with N t atoms can be described 
by the atomic states: 



VN~t- 
Jx _ gt(gt)"*- 1 



Hi) 
Hi) - 



E( = N t ei 
|0) Ei = (JV t - + e 2 (9) 



\^i) = ^== 10) H = ( N t 2) Cl + 2e 2 

where in the ground state Hi), all atoms stay in the low- 
est atomic state \Si). In the first excited state Hi ), onc 
atom is excited to the IS2) state and all the others stay 
in the lowest atomic state. This result is also confirmed 
by the numerical calculations. 

When the collision term can not be neglected, we nu- 
merical solve the Hamiltonian in Eq. (4). In Fig. 6a, the 
calculated energies for the qubit for a large range of £1 
are plotted. The inset of this plot shows the persistent 
currents of qubit states versus fi - The average currents 
(I ac ) of the two qubit states converge to each other as Oo 
increases. 



When the Raman coupling S7o is tuned slowly, the 
qubit state can be manipulated adiabatically. Here 
"slow" means 



|min{£ 2 (r! ) 



£i(^o)}| 2 » 



dO 
~~dT 



(10) 



where dflg/dt is how fast f2 is tuned, and 
mino {_E 2 (^o) — £a(O )} is the smallest energy differ- 
ence between the qubit states during the tuning process. 
As is shown in Fig. 6a, it reaches its minimum at the left 
most end when fi is small. Hence the switching process 
takes a time of milliseconds. 

This adiabatic process can be exploited for efficiently 
initializing the qubit to its ground state. Starting from 
the large limit, we prepare the qubit in its ground 
state \tpi), which is equivalent to preparing all the atoms 
in state \Si) and which can be achieved easily. Then, the 
Raman coupling is adiabatically decreased to the work- 
ing regime so that the ground state |V>g) is reached with 
high fidelity. 



C. Quantum Nondemolition Measurement 



Second, and most important, the adiabatic switching 
provides a scheme for a quantum non-demolition mea- 
surement of the qubit. Starting from the working param- 
eters of the qubit where 2Q N a /Uo = 70 and assuming 
an initial state a^g) + PHD1 ^0 is slowly increased to 
the Rabi regime. When QoN a » Uq, the qubit state 
evolves to a\tpi) + f3\ipi ), a superposition of the states in 
Eq. (9). As the increase of fio is adiabatic, no transition 
to the excited state is induced. Then, a dark-state mea- 
surement scheme is performed on the qubit. Namely, a 
laser pulse is applied that excites the atomic state IS2) 
to an excited state |e) and does nothing to the atoms in 
the states |Si) and IS3). The state |e) emits a photon 
via spontaneous emission which is then detected. As can 
be seen from Eq. (9), when the laser is applied to the 
ground state \if>i), no transition happens and no photon 
is emitted; when the laser is applied to the second state 
\ipi), one atom is excited to the state |e) and one photon 
is emitted. Hence, this approach achieves a projective 
measurement of the qubit. 

The laser pulse applied after the adiabatic switching is 



sU 



(ii) 



where et and e are the operators for the excited state 
|e). It is easy to prove that single atom states S\ and 
S3 are dark states of this operator that can not be ex- 
cited by this pulse (as they are orthogonal states of the 
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Hamiltonian in Eq. (7). We have: H m \ipi) — and 
'KmWl) = \(N - l) Sl , l e )- By performing a single pho- 
ton measurement with the quantum jump approach, the 
probability of the qubit in and hence in \ipf} origi- 

nally, can be detected. For the qubit in its ground state, 
the measurement won't do anything and is a quantum 
nondcmolition measurement. 



VI. DECOHERENCE 

A major obstacle in the pursue of quantum compu- 
tation with solid-state qubits is the strong coupling to 
noise and the resulting low quality factor. In experi- 
ment, the measured decoherence time for the supercon- 
ducting qubits is T 2 = lOOnsec, while the gate time is 
Tgate = lOnsec [8]. 

In the atomic qubit presented in this paper, the quality 
factor due to decoherence is is high compared with that 
of the superconducting qubit. The qubit is designed to 
be insensitive to the major factors that can result in de- 
coherence. For example, all the energies involved in qubit 
operation are much lower than the trapping frequency in 
the longitudinal direction of the trap u>\\ = 3.7 kHz, which 
keeps the atoms in the motional ground state during gate 
operations. Other factors such as the inaccuracy in the 
Raman couplings, the particle loss from the trap and 
the spontaneous emissions can be well neglected within 
a time of seconds. 

The fluctuation of the number of atoms could induce 
severe qubit decoherence when the number of atoms is 
large. For example, the decoherence rate due to single 
particle loss grows linearly with N t and the decoherence 
rate due to three body collision increases with 7V t 3 . Our 
study shows that for single particle loss process with cou- 
pling constant 70, the decoherence rate at N t = 15 is 470, 
and the decoherence due to three body collision can be 
neglected. 

A. Effect of Single Atom Loss 

Consider for example a single atom loss characterized 
by a loss rate 7. The time evolution of the density matrix 
is described by the following master equation 

^ = -i[Hi(i), p% 

<y v-T (12) 

-Qj: = -7o 2J (al^aP* + Alfla ~ 2a a p*fflt ) , 

a 

where p* is the density matrix of the qubit in the inter- 
action picture and the atomic losses in different modes 
are summed up. 

The density matrix can be decomposed into the 
Hilbert spaces of different number of atoms: p* = 
EnPij^NnXinl, where p\f = (i„|p*|j'„) is the element of 



the density matrix with n atoms and \i n ,jn) are qubit 
states of n atoms. Substituting this expression into 
Eq. (12) and assuming an initial density matrix p° with 
Nt atoms, we have 

Pi! f (st) = p^ + p^- 1) 

pW) = p° - St l0 ]T AiA aP ° + P °AlA a (13) 
p (N t -i) =+2Stlo f2A a p°Ai 

a 

where the matrix (Ajj)^- = (iN t \& a \jN t -i)- Starting from 
Nt atoms in the trap, when one atom leaks out, the qubit 
state is a superposition of the eigenstates of (N t — 1) 
atoms. The decoherence rate is slowed down by the fact 
that the remaining system of (N t — 1) atoms largely over- 
laps with the original qubit states in the (N t — l)-atom 
basis. The decoherence rate is expressed as 

7cff = 7 onu«{£<*l4 A a \V) - |(*|4|*)| 2 }. (14) 



Numerical results show that 7 e g- grows linearly with 
number of atoms in the trap as is plotted in Fig. 7a. The 
decoherence is slowed down by a factor of two by the sec- 
ond term in the above equation. At N t = 15, 7 e ^ = 470. 
In the inset of Fig. 7a, the dependence of decoherence on 
the phase of the Raman coupling <fro is plotted which is 
flat in the range of interest. 



B. Three-body Collision Loss 



One of the main decoherence against this qubit is 
the three-body collision loss. The three body process: 
A + A + A — ► A2 + A describes that when three atoms 
collide, two atoms form a bounded molecular state with 
a binding energy of order of h 2 /mag which is several or- 
ders larger than the trapping frequency, where a s is the 
s-wave scattering length. As a result, the molecule and 
atom gain very large kinetic energy after the collision 
and escape from the trap. This process damages the 
coherence of the qubit states. The three-body loss is 
characterized by 7<j 3) = K 3 (3ir 3 )- 3 / 2 p 2 /72N 2 where K 3 
is the three body collision rate in [18] and J d 3 x\tp(x)\ 6 
gives the dependence on density. We apply the same 
approach as that in the single atom loss to calculate 
the effective decoherence rate and the results are plotted 
in Fig. 7b. It is shown that J c s / 7^ N 2 grows linearly 
with N t and at N t = 15, with p = 3 x 10 14 cnr 3 and 
K3 = 10~ 28 cm 6 /sec, we have 7 c ff ~ 10~ 4 which gives a 
very long decoherence time. 
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VII. CONCLUSIONS 
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We have presented a scheme for implementing an 
atomic "flux" qubit with atomic Josephson junctions, 
which are generated by Raman lasers that introduce cou- 
pling between internal modes of atoms. By trapping 
three internal modes and coupling them with the Ra- 
man pulse, a three junction loop is constructed. The 
collision interaction between the atoms provides the ana- 
log of the capacitance energy. With small number of 
atoms, the qubit presents the main features of the meso- 
scopic circuit — superconducting flux qubit: the butterfly 
shaped energy spectrum, the persistent currents and the 
local wave function in phase basis. We have outlined 
methods for the implementation of quantum logic gates 
with fast switching of Raman pulses, the state initializa- 
tion, and we have presented a qubit measurement scheme 
by adiabatic switching of the Josephson coupling and ob- 
servation of quantum jumps. Furthermore, we have given 
detailed analysis of possible imperfection and decoher- 
ence of the qubit. 

The solid-state qubits suffer severely from noise, which 
may become the biggest obstacle in implementing those 
qubits. However, the solid-state proposals are easy to 
scale up and control with existing technology. The qubit 
proposed in this paper inherits many of the merits of the 
superconducting qubits. For one thing, almost all the 
parameters of the qubit can be very well controlled by 
external sources which increases the flexibility of qubit. 
The system is in principle scalable by storing the atomic 
flux qubit in wells of the ID optical lattice. Compared 
with superconducting qubit, the atomic Josephson junc- 
tion qubit has the advantage of not subjecting to severe 
environmental disturbance and having a long decoher- 
ence time. Hence, an array of the atomic qubits can be 
arranged in a space to simulate a "clean" array of su- 
perconducting qubits and perform certain quantum gate 
operations. Clearly, one of the main differences to the su- 
perconducting case is the significantly slower time scale 
of operations. 

In summary, our study shows that the atomic systems 
can be designed to be a clean realization of the Josephson 
junction circuits and keep the merits of exploring macro- 
scopic/mesoscopic degrees of freedom and a long decoher- 
ence time. In this system, the Josephson couplings can 
be controlled with large flexibility by adjusting the power 
and phases of the laser beams. The collision interaction 
can also be adjusted in a large range by magnetic-optical 
means such as tuning around the Feshbach resonances 
[20]. Moreover, the trap geometry and the interaction 
between neighboring qubits can be chosen to suit differ- 
ent experiments. 
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(c) (d) 

FIG. 1. The superconducting flux qubit. (a) the circuit of 
the flux qubit. (b) the potential energy for the qubit. The 
black centers are local maxima and the white centers are local 
minima, (c) the energy and the average current of the qubit 
versus the flux. The arrows indicate the qubit states with 
opposite currents. The double well potentials at the corre- 
sponding flux are plotted, (d) the measurement of the qubit 
by a dc SQUID. 
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FIG. 2. The atomic Josephson junction qubit. (a) Atoms 
trapped in the cigar-shaped optical potential by laser beams, 
(b) Left: the internal modes coupled by Raman pulses. Right: 
the superconducting flux qubit. 
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FIG. 3. The energy and average current of the qubit states 
versus the phase <fr . (a) N t = 15 atoms. Solid lines: for 
symmetric interaction with Uo = 550 Hz; dashed lines: for 
symmetric interaction with Uo = 0; dotted lines: for interac- 
tion given in the text, (b) Energies of the qubit with various 
number of atoms for the symmetric interaction. Solid lines: 
N t = 15; thin dotted lines: N t = 10; dashed lines: N t = 30; 
dotted lines: N t = 50. 
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(a) (b) (c) (d) 

FIG. 4. Contour plots of the probability \ip g ((p a , >fib)\ 2 of 
the ground state wave function of the qubit: (a) Uo = 0.01, 
<f>o = 1/2; (b) Uo = 0.2, O = 1/2; (c) U = 0.01, O = 0.495; 
(d) Uo = 0.2, 0o = 0.495. 
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(c) 

FIG. 5. The absolute value of the elements of the unitary 
transformations for quantum logic gates, \Uij\. The transfor- 
mations are on the lowest six states of the qubits with the 
lowest two states the | |) an d | I) states of the qubit. (a) 
single-qubit NOT gate. The labels indicate the lowest qubit 
states from 1 to 6. (b) single-qubit Hadamard gate. Labels 
are same as in (a), (c) two-qubit gate by a 36-pulse sequence. 
The labels 1 to 4 are the qubit states: | |t, tl, IT, 11)- The 
rest are higher states. 
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(a) (b) 

FIG. 6. Adiabatic switching of Raman tunnelings. (a) The 
energy spectrum of the qubit versus the Raman coupling 0. 
The Raman coupling is plotted in unit of the Raman coupling 
Slo for the designed qubit. The energy differences between 
states AE12 and AE23 are indicated by arrows. The inset 
shows the average current {I}/ SI in the same range of fl. (b) 
The laser pulse H m of the QND measurement after the adi- 
abatic switching. The coupling constants 012, /fe, 72 show the 
relative phase between the three components of the pulse. 
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FIG. 7. Decoherence rate by Eq. (14). (a) Single atom loss 
rate versus number of atoms. Inset: Single atom loss rate 
versus the phase <f>o of Raman coupling Q a b- (b) Three-body 
loss rate in unit of 7oiV t 2 versus the number of atoms. 



